Synchronization on community networks 
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In this Letter, we propose a growing network model that can generate scale-free networks with 
a tunable community strength. The community strength, C, is directly measured by the ratio of 
the number of external edges to internal ones; a smaller C corresponds to a stronger community 
structure. According to the criterion obtained based on the master stability function, we show 
that the synchronizability of a community network is significantly weaker than that of the original 
Barabasi- Albert network. Interestingly, we found an unreported linear relationship between the 
smallest nonzero eigenvalue and the community strength, which can be analytically obtained by 
using the combinatorial matrix theory. Furthermore, we investigated the Kuramoto model and 
found an abnormal region (C < 0.002), in which the network has even worse synchronizability 
than the uncoupled case (C = 0). On the other hand, the community effect will vanish when C 
exceeds 0.1. Between these two extreme regions, a strong community structure will hinder global 
synchronization. 

PACS numbers: 89.75.Hc, 05.45.Xt 



Synchronization is observed in a variety of natural, so- 
cial, physical and biological systems, and has found ap- 
plications in a variety of fields 0. The large number 
of networks of coupled dynamical systems that exhibit 
synchronized states are subjects of great interest. In 
the early stage, these studies are restricted to either the 
regular networks Q, or the random ones Recently, 
inspired by the new discovery of several common char- 
acteristics of real networks, the majority of the studies 
about network synchronization focus on networks with 
complex topologies. The effects of average distance 
heterogeneity [5|, clustering [6(, and weight distribution 
on network synchronizability have been extensively 
investigated. 

Besides the small-world and scale-free properties, it 
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FIG. 1: (Color online) The inverse of the smallest nonzero 
eigenvalue 1/Ai (left) and the eigenratio R (right) vs C. The 
red circles, black squares, and green triangles represent the 
cases of n = 2, n = 5, and n = 10, respectively. As shown in 
the insets, the same data can be well fitted by a straight line 
with slope ~ — 1 in the log- log plot, indicating the relation 
j- ~ C~ x . All the data are obtained as the average over 10 
realizations, and for each realization, the network parameters 
N = 5000 and m — 3 are fixed. 



has been demonstrated that many real networks have 
the so-called community structure 0. Qualitatively, a 
community is defined as a subset of nodes within a net- 
work such that connections between the nodes therein 
are denser than that with the rest of the network 
Very recently, by applying the epidemiological models on 
community networks, it was found that the network epi- 
demic dynamics are highly affected by the community 
structure [T^j- To date, however, the issue of synchro- 
nization on community networks has not been fully in- 
vestigated. Based on a toy network model with a tunable 
community strength, in this Letter we intend to provide 
a first analysis on how community structure affects the 
network synchronizability. 

Our model starts from n community cores, each core 
contains mo fully connected nodes. Initially, there are no 
connections among different community cores. At each 
time step, for each community core, one node is added. 
Thus, there are in total n new nodes being added in one 
time step. Each node will attach m edges to existing 
nodes within the same community core, and simultane- 
ously ml edges to existing nodes outside this community 
core. The former are internal edges, and the latter are 
external edges. Note that, the m and ml are not neces- 
sary to be integers, for example, to generate 2.7 edges can 
be implemented as follows: Firstly, generate 2 edges, and 
then generate the third one with probability 0.7. Similar 
to the evolutionary mechanism of Barabasi- Albert (BA) 
networks we assume the probability of choosing an 
existing node i to connect to is proportional to z's degree 
hi. Each community core will finally become a single 
community of size N c , and the network size N = nN c . 
By using the rate-equation approach |l2f. one can eas- 
ily obtain the degree distribution of the whole network, 
p(k) r~j fc -3 . For simplicity, we directly use the ratio of 
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external edges to internal ones, C = to measure the 
strength of the community structure. Clearly, a smaller 
C corresponds to sparser external edges thus a stronger 
community structure. 

Next, we investigate how the community strength C 
affects the network synchronizability. Consider N identi- 
cal dynamical systems (oscillators) with the same output 
function, which are the nodes of a network and coupled 
linearly and symmetrically with neighbors through edges. 
The coupling fashion ensures the synchronization mani- 
fold be an invariant manifold, and the dynamics can be 
locally linearized near the synchronous state. The state 
of the ith oscillator is described by x l , and the set of 
equations of motion governing the dynamics of the iV 
coupled oscillators is 

JV 

x^FfxV^GyH^), (1) 
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where x 4 = F(x J ) governs the dynamics of the ith indi- 
vidual oscillator, H(x. J ) is the output function, a is the 
coupling strength, and the N x N Laplacian G is given 
by 

!ki for i = j 
-1 forjeA, (2) 
otherwise. 

Being positive semidefinite, all the eigenvalues of G are 
nonnegative reals and the smallest eigenvalue Ao is always 
a single zero, for the rows of G have zero sums. Thus, the 
eigenvalues can be ranked as 0=Ao < Ai < • • • < Ajv-i- 
If the synchronized region is left-unbounded, according 
to the Wang-Chen (WC) criterion the network syn- 
chronizability can be measured by the inverse of the 
smallest nonzero eigenvalue 1/Ai: the smaller the 1/Ai, 
the better synchronizability, and vice versa. We show 
the numerical results about the relation between 1/Ai 
and C in the left plot of Fig. 1. Clearly, the community 
structure will hinder the global synchronization. More 
interestingly and significantly, we found an unreported 
linear relationship, Ai cx C . This seems a universal law 
for community networks if the community structure is 
sufficiently strong (i.e. C is sufficiently small). 

Denote by Gi the N c x N c Laplacian of the ith commu- 
nity, Hi the submatrix of G consisting of all the rows and 
columns corresponding to the nodes of the ith commu- 
nity. Note that Hi and Gi are different only in diagonal 
elements, and H^s smallest eigenvalue A^ is positive. 
By using the combinatorial matrix theory, one can first 
prove that the smallest nonzero eigenvalue Ai of G is 
equal to the minimal one of all the smallest eigenvalues 
\f Q (i = 0, 1, • • • ,n): Ai = mini<j<„A^. And, then, one 
is able to prove that for each matrix Hi, A^ is approxi- 
mately linearly correlated with the community strength 
C. The strict and full proof is fairly complicated and is 
omitted here. 

If the synchronized region is finite, according to the 
Pecora-Carroll-Barahona (PCB) criterion the net- 
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FIG. 2: (Color online) Order parameter M vs coupling 
strength a for different values of the community strength C. 
The solid line represents M = 0.403. All the data are obtained 
as the average over 100 realizations. For each realization, the 
network parameters N = 5000, n = 5, N c — 1000 and m = 3 
are fixed. 




FIG. 3: (Color online) Order parameter M vs community 
strength C for different values of the coupling strength a. The 
dash line represents M = 0.403. All the data are obtained by 
the average over 100 realizations. For each realization, the 
network parameters N = 5000, n — 5, N c — 1000 and m = 3 
are fixed. 



work synchronizability can be measured by the eigenratio 
R = AjV-i/Ai: the smaller it is, the better synchroniz- 
ability will be, and vice versa. We also have checked that 
the maximal eigenvalue is not sensitive to the change of 
the community strength (it slowly diminishes as C de- 
creases), thus both the WC and PCB criteria give qual- 
itatively the same result (see the right plot of Fig. 1 for 
the case of the PCB criterion). 
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Hereinafter, we investigate a network of nonidentical 
Kuramoto oscillators 0], obeying the coupled differen- 
tial equations 



d(p t 
~dt 



(3) 



where < <f>i < 2n are phase variables, u>i are intrinsic 
frequencies, i = 1, 2, • • • , N, a is the coupling strength, 
and [aij] is the adjacency matrix (aij = 1 iff nodes i and j 
are connected) . Initially, (pi and are randomly and uni- 
formly distributed in the intervals [0,2-71") and [—0.5,0.5], 
respectively. The numerical results are obtained by in- 
tegrating Eqs. (3) using the Runge-Kutta method with 
step size 0.01. To characterize the synchronized states, 
we use the order parameter 
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where {•} signifies the time averaging. The order pa- 
rameters are averaged over 10 4 time steps, excluding the 
former 5000 time steps, to allow for relaxation to a steady 
state. Clearly, M is of order 1/V~N if the oscillators are 
completely uncorrelated, and will approach 1 if they are 
all in the same phase. 

Fig. 2 reports the simulation results for different com- 
munity strength C. Whatever the value of C, the param- 
eter M increases sharply after a critical point a c = 0.6. 
This point is just the same as the critical point at which 
synchronized behavior emerges in each separate commu- 
nity (BA network of size 1000 and with average degree 
6). For all the cases with C > 0.003, a strong community 
structure (i.e. a smaller C) hinders global synchroniza- 
tion. It is difficult to harmonize different communities 
based only on a very few external edges. The community 
effect becomes lower as C increases. For sufficiently large 
C (see the cases of C = 0.15 and C = 0.30 in Fig. 2), 
the network synchronized behavior is almost the same 
as that of the original BA network; that is, community 
effect vanishes. 

Consider a network consisting of n uncoupled commu- 
nities (C = 0), and each community itself can approach 
a nearly completely synchronized state. The order pa- 
rameter M of the whole network is 



M(n) 
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n(2n) r 
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where x = \l Ei si 11 ^) 2 + Ei cos</>i) 2 . The numerical 
result M(5) 0.403 is represented by a horizontal line in 
Fig. 2. The corresponding simulation with C = shows 
that M slightly fluctuates around 0.403 for sufficiently 
large a, which is in accordance with the above numer- 
ical result. Interestingly, we found an abnormal region 
C < 0.002, in which the networks have even worse syn- 
chronizability than the uncoupled (unconnected) case. 
Within this region, each separate community can not get 
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FIG. 4: (Color online) The order parameters vs time for the 
whole network and each community. The network parameters 
are N = 5000, n = 5, N c = 1000 and m = 3. The coupling 
strength a — 5 is fixed. 



harmonized with other communities through its few ex- 
ternal edges. On the contrary, the input signals contain- 
ing by these edges disturb the synchronizing process of 
this community. 

Fig. 3 shows the order parameter M vs C for different 
a. As mentioned above, when C < 0.002, M can not 
reach the dash line (M=0.403). We have also checked 
that even for very large tr, M is always smaller than 0.403 
if C < 0.002. The distinct difference between the orig- 
inal BA networks and the present community networks 
vanishes when the density of external edges exceeds 0.1. 
The division of those three regions is qualitative, and the 
borderlines between neighboring regions can not be ex- 
actly determined. However, it provides a clearer picture 
about the effect of the community structure. 

To further understand the underlying mechanism of 
synchronization on community networks, we investigate 
the partial synchronization within a separate community. 
For the ith community, the corresponding order param- 
eter ML is defined as 



(6) 



where the sum goes over all the nodes belonging to the 
ith community. Fig. 4 exhibits the temporal behaviors 
of order parameters for the whole network and for each 
community. The four plots correspond to the cases of 
G = 0.001, C = 0.003, C = 0.02, and C = 0.15, re- 
spectively. In the abnormal region (C = 0.001), due to 
the external disturbance, the order parameter of each 
community is remarkably below 1 even in the long time 
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limit. After Mj reaches its steady value (t w 4), M 
also becomes steady, indicating that the external edges 
can not harmonize different communities, but only intro- 
duce some noise thus hinder the partial synchronization. 
When C > 0.003, a separate community can approach a 
nearly completely synchronized state, and can harmonize 
with other communities; therefore, M will continuously 
increase after Mj gets steady. For sufficiently large C 
(see the case of C = 0.15), the whole network can ap- 
proach the nearly completely synchronized state almost 
as quickly as the single separate community, and the ef- 
fect of the community structure on network dynamics (at 
least for the Kuramoto model) is hardly observed. 

In conclusion, we have proposed propose a scale-free 
network model with a tunable community strength and 
studied its synchronization phenomenon. We found an 
unreported linear relationship Ai oc C, which is the first 
quantitative formula that describes the synchronizabil- 
ity of community networks. We have also checked that 
the maximal eigenvalue is not sensitive to the change 
of the community strength, thus both the Wang-Chen 
and Pecora-Carroll-Barahona criteria give qualitatively 
the same result: The stronger the community struc- 
ture, the worse the synchronizability. Furthermore, we 
have investigated the Kuramoto model in community 
networks. Interestingly, we found an abnormal region, 
in which the networks have even worse synchronizability 
than the uncoupled case. Due to the complicacy of the 
scale-free structure itself, we are unable to give a the- 
oretical and analytical explanation about this observed 
phenomenon. An approximate analytic solution of a sim- 



ilar phenomenon may be obtained based on a more ideal 
community network model, where each community is a 
complete graph (Very recently, this ideal model has also 
been investigated to show the effect of modular number 
on network synchronization |l6l]h Beyond this abnor- 
mal region, analogous to the result from the approach 
of the master stability function, increasing the density of 
external edges will sharply enhance the network synchro- 
nizability. However, when the density of external edges 
exceeds 0.1, the synchronized behavior becomes almost 
the same as that of the original BA networks and fur- 
ther enhancement can not be achieved. This result is not 
only of theoretical interest, but also significant in prac- 
tice if one wants to enhance the synchronizability of com- 
munity networks by adding external edges. Finally, we 
would like to point out that, although the present model 
is very simple, it provides a useful way to have detailed 
understanding about the effect of the community struc- 
ture on network dynamics since the community strength 
in the model is adjustable. We believe that this model 
can also be applied to the studies on many other network 
dynamical processes. 
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